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ABSTRACT 

, Among dozens young pulsar wind nebulae, some have been detected in TeV 7-rays 

(TeV PWNe), while others have not (non-TeV PWNe). The TeV emission detectability 
is not correlated either with the spin-down power or with the characteristic age of their 
^ c| central pulsars, and it is an open problem what determines the detectability. To study 

Q_f this problem, we investigate spectral evolution of five young non-TeV PWNe, 3C58, 

' ; G310.6-1.6, G292.0+1.8, Gll.2-0.3 and SNR B0540-69.3. We use a spectral evolution 

^ ■ model which has been developed to be applied to young TeV PWNe in our previous 

works. TeV 7-ray flux upper limits of non-TeV PWNe give upper or lower limits on 
\ parameters, such as the age of the PWN and the fraction of the spin-down power 

going to the magnetic energy injection (the fraction parameter) . Combined with other 
independent observational and theoretical studies, we can guess a plausible value of 
the parameters for each object. For 3C58, we prefer the parameters with an age of 
2.5 kyr old and the fraction parameter of 3.0 x 10~ 3 , although the spectral modeling 
I alone does not rule out a shorter age and a higher fraction parameter. The fraction 

. parameter of 3.0 x 10 -3 is also consistent for other non-TeV PWNe and then the value 

f — 1 is regarded as common to young PWNe including TeV PWNe. Moreover, we find that 

the intrinsic properties of the central pulsars are similar, 10 48_50 erg for the initial 
rotational energy and 10 42_44 erg for the magnetic energy (2 x 10 12 — 3 x 10 13 G for 
the dipole magnetic field strength at their surfaces). The TeV detectability is correlated 
with the total injected energy and the energy density of the interstellar radiation field 
around PWNe. Except for G292.0+1.8, a broken power-law injection of the particles 
well reproduces the broadband emission from non-TeV PWNe. 
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1 INTRODUCTION 

Rotation-powered pulsars release most of their rotational 
energy as a relativistic outflow of a magnetized electron- 
positron plasma called pulsar winds. The pulsar wind col- 
lides with a surrounding supernova ( SN) ejecta and cre- 
ates a pulsar wind nebu la (PWN) (|Rees fc Gunnl 1 1974 
iKennel fc Coronitil Il984ah . The resultant PWN is a cloud 
of the relativistic magnetized electron-positron plasma dis- 
sipated at a termination shock of the pulsar wind and shines 
from radio through TeV 7-rays via synchrotron radiatio n 
and inverse Compton scattering |Kennel fc Coronitilll984bT ). 

More than fifty PWNe have been identified till now. 
Although they share some common observational character- 
istics, they also show individualities, such as an observed 
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size and a TeV 7-ray to X-ray flux ratio. The differences 
may be due to the properties of the central pulsar (e.g., the 
spin-down power and the surface magnetic field), of the en- 
vironment (e.g., SN ejecta, the interstellar medium and the 
interstella r radiation field (ISRF)) and of the age of each ob- 
ject (e.g.. iTanaka fc Takahara! l20ld. [2011 hereafter TTlil ; 
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iBucciantini et al.ll201ll : lGelfand et al.ll2009t ). In other words 
we can infer the properties of the central pulsars, the en- 
vironment of the objects and the age of the objects from 
observed characteristics. 

The goal of this paper is to understand which prop- 
erties of PWNe affect TeV 7-ray properties. The ground- 
based 7-ray telescopes (e.g., HESS, MAGIC, VERITAS 
etc.) have detected TeV 7-rays from many PWNe (TeV 
PWNe), but some hav e not been detected (non-TeV PWNe). 
iMattana et al.l |2009l ) studied correlations between the TeV 
7-ray luminosity from PWNe and their pulsar properties 
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(the spin-down power and the characteristic age). However, 
they found no correlation (c.f., see their Figure 1). Therefore, 
we expect that other properties characterise the detectabil- 
ity of TeV 7-rays from young PWNe. In this paper, we study 
this problem from the broadband spectrum of PWNe. 

We built a spectral evolution model of young TeV 
PWNe in our previous works and found that they share 
common characteristics, such as a magnetic energy fraction 
of the total energy injected from the central pulsar (TT11). 
Here, we apply the same spectral evolution model to young 
non-TeV PWNe and study their differences and similarities. 
After the applications of our model to non-TeV PWNe, we 
discuss the properties of the whole class of young PWNe. 
Note that our spectral evolution model is applicable only 
to young PWNe which show no signature of an interac- 
tion with a reverse shock of a supernova remnant (SNR). 
This is because our model does not consider complex dy- 
namical behaviors of PWNe in the inte raction phase (e.g., 
iBlondin et al. I l200ll ; iGelfand et alj|2009h . 

We apply the model to five young non-TeV PWNe 
3C58, G310.6-1.6, G292. 0+1.8, Gil. 2-0.3 and SNR B0540- 
69.3 (hereafter B0540-69.3). We pick up these five PWNe 
according to four criteria. (1) They reveal non-thermal spec- 
tra at least in radio and X-rays. (2) They have a flux upper 
limit in TeV 7-rays. (3) They have a central pulsar with a 
known period and its derivative. (4) They have an almost 
spherical shape with a known angular extent, i.e., they are 
young enough not to reveal any signatures of an interaction 
with the SNR reverse shock. 

In Section [2j we describe the fitting procedure which is 
different from that for TeV PWNe. The spectral evolution 
model is the same as our previous works and is summarized 
in Appendix IA1 In Sections [3] and [4] we apply the model to 
the five non-TeV PWNe and present results. Comparative 
discussions about young non-TeV PWNe and about young 
TeV PWNe are made in Section [5] Conclusions are made in 
Section IS] 



2 FITTING PROCEDURE 

For each object, we adopt the distance to the object d, the 
angular extent of the object, and the central pulsar param- 
eters, the period P, its derivative P and the braking index 
n from observations. Then, we can determine the radius of 
the object Rpwn and the energy densities of the ISRF Um 
and Uopt- We have already studied the dependence of our 
calculations on d, Uib, and [/opt (see Section 2.2 of TT11), 
and then we do not consider their uncertainties in this pa- 
per. Adopted parameters for each PWN are listed in Table[T] 
where we also show results obtained in our previous papers 
(TeV PWNe) for comparison. 

Fitting parameters are the age t ag e, the fraction param- 
eter 77 (riL S pi n (t) is the magnetic power injected from the 
pulsar), and the parameters of the particle distribution at 
injection Qi n j(7,i) in Equation (| A2p . 7 max , 7b, 7min, Pi and 
P2. Four of seven parameters, 7 max , 7min, Pi and p2, de- 
termine a broadband spectral shape, and they are usually 
fitted after the determination of the other three parameters. 
Note that only a lower limit of 7 max and an upper limit of 
7min can be deduced from observations in radio and X-rays 
(see Section 2.3 of TT11). It should also be noted that we 



do not fit the observed spectral index in X-rays, but fit the 
flux value. Because spatial variabilities of the X-ray spectral 
index are observed for many PWNe, the one-zone spectral 
model needs not to rep roduce the spectral index in X-ray 
jBucciantini et alj|201ll . TTU). 

Three parameters i agc , f] and 7b are fitted as follows, 
bearing in mind that typical spectra of young PWNe re- 
quire pi < 2 < P2 and 77 <C 1. Basically, for young TeV 
PWNe, three observed values, the power of TeV 7-rays 
-Pic/isrf <x 7b-Etot(tage), the power of radio and X-rays 
fsyn oc 7b??-£tot(iage) and the spectral break frequency ia> oc 
7b?7 1,/2 -B t J / t 2 (f a ge) (note that this is not the synchrotron cool- 
ing break frequency), can be used to fit t ago (or -Etot(tage)), 
77 and 7b. However, for non-TeV PWNe, we have only an 
upper limit in TeV 7-ray, and then we will get upper limits 
for t age and 7b and a lower limit for 77, i.e., one of these three 
parameters should be chosen in other ways. We compare the 
obtained range of t age with other independent studies about 
t agc (see below). Note that i5 t ot(£ age ) is an increasing func- 
tion of t age and then an uncertainty of i agc propagates to 77 
and 7b as 77 oc _E t ~t /3 (t age ) or 7 b cc E^(t^ gc ) to retain the 
observed values of P syn and Vh ■ 

There are several methods to study the age t ago of 
PWNe. Observationally, the age t ago is estimated from an 
expansion rate of PWNe assuming the rate being constant, 
from a Doppler shift of emission lines also assuming the ex- 
pansion velocity being constant, from a proper motion of 
the central pulsar, and from a possible coincidence of his- 
torical events of SNe. Strictly speaking, although the line 
emission is from a SN ejecta component, but not from a 
PWN, this velocity is consistent with the expansi on velocity 
of the PWN for the Crab Nebula. Theoretically, IChevalier! 
(|2005. hereafter C05l ) studied dynamical evolution (expan- 
sion) of PWNe inside a SNR and estimated the age of some 
PWNe taking account of the spin-down evolution of their 
central pulsars and of the SN-types, i.e., the density pro- 
file of the SN ejecta. Because these estimates are not always 
consistent with each other, we just refer to these studies and 
discuss each object individually. 

Once the seven parameters are fitted, then we obtain 
the expansion velocity vpwn, the magnetic field strength 
Bnow and the parameters of a central pulsar, the spin-down 
time ro and the initial rotational energy Lotq. Note that we 
obtain upwn from i age and the observed extent of a PWN, 
i.e., our dpwn is independent of the study by C05 described 
above. We also get the pair multiplicity inside a pulsar mag- 
netosphere k and the bulk Lorentz factor of a pulsar wind 
r w assuming that the pulsar wind is a cold electron-positron 
plasma flow. Fitted, and derived parameters for each PWN 
are also listed in Table Q] Note that an uncertainty of t age 
also propagates to these six derived parameters in the way 
they decrease with increasing t ago except for Lqtq. Depen- 
dence on t age is easy to find from comparing two cases (tage = 
2.5kyr and lkyr) of 3C58 in Table [T] for example. 



3 3C58 

3C58 is a well known filled center SNR, and non-thermal ra- 
diation is observed in radio, infrared and X-rays. Although 
3C58 has been observed in GeV and TeV 7-rays, only flux 
upper limits are obtained. The angular extent is ~ 6' x 9' 
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Table 1. Adopted, Fitted, and Derived Parameters of Young PWNe 



Symbol 


TeV PWNe" 


3C58 


3C58 6 


G310.6-1.6 


G292.0+1.8 


Gil. 2-0.3 


B0540-69.3 


Model 




1 


2 










Adopted Parameters 


a!(kpc) 


2.0 - 13 


2.0 


2.0 


7.0 


6.0 


5.0 


50 


fipwN,now(pc) 


1.0 - 3.8 


2.0 


2.0 


1.3 


3.5 


0.97 


1.2 


P(msec) 


33.1 - 136 


65.7 


65.7 


31.2 


135 


65.0 


50.7 


P(lCr 13 sec ■ sec" 1 ) 


1.56 - 7.51 


1.94 


1.94 


0.389 


7.47 


0.440 


4.79 


n 


2.51 - 3.0 


3.0 


3.0 


3.0 


3.0 


3.0 


2.1 


(7l R (eVcm- 3 ) 


1.0 - 2.0 


0.3 


0.3 


0.3 


0.3 


1.2 


0.5 


[/0PT(eVcm~ 3 ) 


0.5 - 15 


0.3 


0.3 


0.3 


0.3 


2.0 


0.5 


Fitted Parameters 


?f(i(r 3 ) 


1.0 - 8.0 


3.0 


60 


3.0 d 


3.0 d 


3.0 d 


3.0 d 


i age (kyr) 


0.95 - 4.5 


2.5 


1.0 


0.6 


2.7 


2.0 


0.7 


7S la x(10 9 ) 


0.8 - 7.0 


1.0 


0.5 


1.0 


1.5 


1.0 


0.8 


7b (10 5 ) 


0.4 - 6.0 


0.9 


0.5 


30 


0.3 


1.0 


100 


7^ n (10 3 ) 


0.1 - 20 


0,1 


0.1 


5.0 




2.5 


4.5 


Pi 


1.0 - 1.5 


1.0 


1.0 


1.4 




1,5 


1.7 


P2 


2.5 - 2.6 


3.0 


3.0 


3.0 


2.5 


2.6 


2.5 


Derived Parameters 


«PWN( knl sec 


770 - 1800 


780 


2000 


2120 


1270 


470 


1680 


S now (mG) 


6.7 - 85 


17 


40 


17 


16 


17 


40 


fo(kyr) 


0.7 - 3.9 


2.9 


1,1 


11.5 


0.16 


21 


2,1 


L r (10 48 erg) 


2.6 - 74 


8.6 


5.6 


21 


19 


5.2 


22 


ft c (10 4 ) 


3.4 - 420 


23 


15 


3.8 


5.3 


6.3 


9.5 


r w (io 5 ) 


0.07 - 2.7 


0.31 


1.6 


1.8 


0.90 


1.0 


1.8 



a Taken from Table 1 of TT11 excluding the models whose r\ are significantly different from a 
few X 10~ 3 (model 1 of G21.5-0.9 and G54.1+0.3 and both models of Kes 75). 
^ The case considering SN 1181 (t age ~ Ikyr). 



The values are upper limits for 7 m i n and T w 
We fit the broadband spectrum assuming ?/ 



and lower limits for 7 max and / 
= 3.0 X 10~ 3 for these objects. 



in radio (|Bietenholz et all l200ll ; iBietenhohd I2006J ). A cen- 
tral pulsar of 3C58 (PSR J0205+6449) is observed in ra- 
dio, X-rays, and 7-rays with a period P — 6.57 x 10 _2 s, 
and a time derivativ e P = 1.94 x 10~ 13 s s" 1 (r c = 5.4kyr) 
dCamilo et al.ll2002bl ; iMurrav et ai1l2002l ; iLivingstone et all 
I2009T ). We assume n = 3. The age has bee n estimated with 
three different methods. iFesen et all (|20Qgf ) showed an aver- 
age expansion velocity upwn ~ 770km sec -1 measured from 
infrared emission lines. Corresponding a ge is t agG ~ 2.5kyr 
for -Rpwn ~ 2pc. IStephenson fc Greenl <|2002h argued that 
3C58 might be a historical SN in 1181 A.D., i.e., i agc ~ lkyr. 
An expansio n rate of a radio synchrotron nebula suggests 
i a gc ~ 7kyr (|Bietenhold 12006 ]) but i agc > r c is theoretically 
untenable and cannot be treated in our model because the 
spin-down time becomes to < 0. Moreover, the correspond- 
ing expansion velocity is too small as a young PWN. We thus 
study th e first two case s. We adopt that the distance to 3C58 
is 2 kpc (|Kothesll20T0l ) and then we approximate 3C58 as a 
sphere of radius i?pwN ~ 2pc. For the ISRF energy densi- 
ties, we adopt (i/m, [/opt) = (0.3eVcm -3 , 0.3eVcm -3 ). 




Figure 1. Model spectrum of 3C58 for t ago = 2.5kyr. 
The observational data, the upper limits and the sensitiv- 



ity of CTA (50 hours) are also plotted. The observed data 


are taken from 


IWeiland et all J201lll ISalter et al.l Jl98Sfl: 


lAMI Consortium 


: Hurlev-Waker et all ll2009f);lKothed l|2006h (ra- 


dio), Slane et al. 


feOOSf) (infrared), Torii et alj J200of) (X-rav). 


lAckermann et al 


ll201lll lAliu et al.1 J2008h; lAnderhub et al.l 



Used parameters are tabulated in Table [T] 
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Figure 2. Model spectrum of 3C58 for t agc = l.Okyr. The ob- 
servational data are the same as Figure [T] Used parameters are 
tabulated in Table [T] 



3.1 Results 

First, we show a model spectrum of t agc = 2.5kyr (Figure [T]). 
We reproduce the observed non-thermal spectrum with pa- 
rameters 77 = 3.0 x 10~ 3 , 7 max = 1.0 x 10°, -y b = 9.0 x 10 4 , 
7min = 4.0 x 10 3 , pi = 1.0, and p 2 = 3.0. IC/IR and 
IC/CMB dominate in 7-rays. Because the predicted TeV 
7-ray flux in Figure[T]is almost comparable with the observa- 
tional upper limit, we can put an upper limit of t agc ^ 2.5kyr 
and a lower limit of 77 ^ 3.0 x 10 -3 . Accordingly, derived 
parameters are LqTq = 8.6 x 10 48 erg, wpwn = 780km sec -1 , 
B now ^ 17/jG and to ^ 2.9kyr. The fitted, and derived pa- 
rameters are tabulated in Table [T] 

Figure [2] shows a model spectrum with another param- 
eter set being able to reproduce the observations. For this 
case, supposing that 3C58 is associated with SN 1181, we 
take tagc = l.Okyr. The predicted 7-ray flux in Figure [2] is 
significantly smaller than that in Figure [T] From Table [I] 
we can see how the uncertainty of t ag e propagates to the 
fitted and derived parameters. We find the fraction param- 
eter being an order of magnitude larger 77 = 6.0 x 10~ 2 
(Bnow = 40/xG) and the total injected energy is five times 
smaller E to t = 1.1 x 10 48 erg than those when t agc is taken 
to be 2.5 kyr. Note that here the fraction parameter does 
not change as 77 oc E to ^ 3 (t ag e) as was discussed in Section [2] 
because adiabatic and synchrotron coolings of particles are 
significant for this short age. 



tagc = 2.5 kyr (Figurefl]). Future observations in TeV 7-rays 
will distinguish the models. 

Although past dynamical and spectral evolution stud- 
ies of 3C58 also obtained a s imilar age of 2.4kyr (C05) and 
2kyr (|Bucciantini et al]|201lf ). they adopt different distance 
d = 3 .2kpc while we adopt d = 2kpc taken from iKothesI 
|2010] ). Despite thi s difference, we try to co mpare our model 
with the study by iBucciantini et al. (|201ll l. They obtained 



tagc ~ 2kyr and B(f agc ) ~ 43/^G. Because the magnetic 
field strength is determined from the observed ratio of syn- 
chrotron and inverse Compton scattering Psyn/-Pic/iSRF> it 
does not depend on the distance to 3C58. Their magnetic 
field strength ~ 43/jG is consistent within our condition 
S now ^ 17/kG and -Pic/isrf is similar to Figure [2] (-Bnow = 
40/jG). 

Radio observat ions of 3C58 suggested a flux evolution. 
Ilvanov et al.l (|2004l ) found non- stationary variations of the 
flux densities at 1.5 — 8.5 GHz. Although they found a fre- 
quency independent increase of the radio flux ~ 0.3%yr~ 
during the period 1965 — 1986, they also found a frequency 
dependent variation during the period 1966 — 2003, where 
both increase and decrease of fluxes are found ranging from 
-0.3 to 0.5 %yr _1 . Such a cha otic behavior is different from 
the case of the Crab Nebula |Vinvaikinll2"007l ) and is unex- 
pected as a PWN. On the other hand, a recent observed rate 
of variabilities bv lWeiland et alj (|201ll ) (seven -year WMAP 
observations) seems to be zero within the margin of error, 
-0.05±0.09%yr _1 at 23GHz, 0.14 ± 0.15%yr _1 at 33GHz, 
-O.^iO.lS^oyr" 1 at 41GHz, -0.43±0.35%yr _1 at 61GHz, 
and 0.47 ± 0.80%yr _1 at 94GHz. In our model (t agc = 
2.5kyr), the current decrease rates are about — 0.04%yr~ 
at 1GHz, -0.05%yr _1 at 10GHz, -0.06%yr _1 at 100GHz. 
They are almost independent of frequencies and are a factor 
of four smaller than the case of the Crab Nebula. Because 
the absolute radio flux of 3C58 is also two orders of mag- 
nitude less than that of the Crab Nebula, we conclude that 
such a small rate of decrease is difficult to detect and that 
our result is consistent with WMAP observations. 

The high energy power-law index at injection p 2 = 3.0 
for 3C58 is different from young TeV PWNe, for which typ- 
ically p 2 ~ 2.5 (TT11). When we use p 2 ~ 2.5, the cal- 
culated X-ray flux becomes almost an order of magnitud e 
larger than the observed X- ray flux by iTorii et alj (|2000T ). 
This result does not change even if we ignore the Spitzer 
observation in infrared. 



3.2 Discussion 

Both spectra in Figures [T] and [2] reproduce the observations 
from radio through X-ray and are compatible with the 7- 
ray upper limits. Although an assoc iation with SN 11 81 is 
suggested, some observational (e.g., iFesen et al.ll200si ) and 
theoretical (e.g., C05) studies discuss t ag c > 2.0kyr. In ad- 
dition to these studies, it is reasonable that the fraction 
parameter 77 = 3.0 x 10~ 3 is similar to young TeV PWNe, 
typically 77 = a few x 10 -3 (TT11). We conclude that 3C58 
is about 2.5 kyr old. For example, the fraction parameter 
becomes 77 ~ 10 -2 for t agc = 2.0 kyr. The predicted TeV 
flux for tagc = 2.0 kyr is about twice as small as that for 



4 OTHER YOUNG NON-TEV PWNE 

We also study other four young non- TeV PWNe, G310.6-1.6, 
G292.0+1.8, Gil. 2-0. 3 and B0540-69.3. We found that their 
flux upper limits in TeV 7-rays only weakly constrain the 
fitting parameters. For example, we obtain a lower limit of 
i)^5x 10 -4 for G310.6-1.6 and lower limits of 77 are much 
lower than 10 for other three objects. Hence, 77 being a 
few x 10 -3 is also consistent for these four PWNe as well 
as 3C58 and young TeV PWNe. We regard that 77 ~ a few 
x 10~ 3 is reasonable for all of young PWNe. Figures [3] — 
©show model spectra of G310.6-1.6, G292.0+1.8, Gil. 2-0.3 
and B0540-69.3 when the fraction parameters are taken to 
be 77 = 3.0 x 10~ 3 together with observational data. Adopted, 
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Figure 3. Model spectrum of G310.6-1.6. The observational 
data, the upper limits and the sensitivity of C TA (50 hours) 
are al so plotted. The observed data are taken from lRenaud et al,l 
||2010| ) (radio, X-ray and 7-ray). Used parameters are tabulated 
in Table [T] 



fitted, and derived values corresponding to Figures [3] — [S] are 
listed in Table[T] Below, we briefly discuss each object. 



4.1 G310.6-1.6 

G310.6-1.6 is a recently discovered composite SNR at d = 
7kpc, and the associated pulsar (PSR J1400-6325) has 
= 5.1 x 10 3 7 erg s _1 and r c = 12.7kyr (n = 3) 



(|Renaud et al.l2010 ] ). A model spectrum with observations is 
shown in Figure(3] Renaud et al.l (|2010l ) discussed a one-zone 
leptonic model of G310.6-1.6, and they obtained conditions 
of the magnetic field (> 6/iG) and the age (< 1.9kyr) from a 
TeV 7-ray flux upper limit. Our results are compatible with 
their conditions. 

For 77 = 3.0 x 1(T 3 (B now = 6.7/iG), the age becomes 
f aE c = 0.6 kyr and then G310.6-1.6 is one of the youngest 
We consider that a larger r],i.e., a larger upwn, is 



''age — 

PWNe 
unlikely because wpwn 



2120km sec" 1 for t a 



0.6 kyr is 



already large compa red with other young PWNe. Moreover, 
iRenaud et ail |2010l ) suggested t agc < lkyr based on some 
theoretical studies of dynamical evolutions of a composite 
SNR. Our result of t agc = 0.6 kyr shows that the predicted 
TeV 7-rays from G310.6-1.6 is about three times smaller 
than the HESS upper limit. 



4.2 G292.0+1.8 

G292.0+1.8 is a well-known Cas A-like (O xygen-rich) SNR 
(IPark et al. l l2007l ; lGhavamian et al]|2012t) and has the cen- 
tral pulsar (PSR J1124-591 6) with L spin = 1.1 x 10 37 erg s" 1 
and r c = 2 9kyr (n = 3) dCamilo et"al.|[2002al ). We adopt 
d — 6kpc jGaensler fc Wallace! 20031 ) . A model spectrum 
with observations is Figure [4] Although optical and infrared 
observations found a compact torus-like feature aroun d PSR 
J1124-5916 (|Zharikov et alj|2008l ; IZvuzin et al.ll2009t ), these 
do not seem to account for the total PWN emission at these 
wavelengths (plotted as a lower limit in Figure 2]). This is 
because the detected structure is much smaller than radio 
and X-ray nebula. 
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Figure 4. Model spectrum of G292.0+1.8. The observa- 
tional data, the upper limits and the sensitivity of CTA (50 
hours) are also plo t ted. The obser ved data are taken from 
iGaensler fc Wallace! ll2003l) (radio). IZvuzin et all l|2009l) (in- 
frared). IZharikov et alj <2008ft (op tical). iHuehes et alJfeOOll) (X- 
ray) and lAckerrn^imetal (7-ray). Used parameters are 

tabulated in Table [T] 



No other broadband spectral study of G292.0+1.8 is 
found in the literature to o ur knowledge. Althou gh C05 ob- 
tained t ago ~ 3.2kyr and IWinker et al.l (|2009h estimated 
tagc ~ 3.0kyr from a proper motion of oxygen-filaments, 
tagc > t c is forbidden with n = 3 (C05 considered the 
case n < 3). Our result is t ag o = 2.7kyr < r c . Although 
the values of t age and n are somewhat different between 
our and C05's models, the internal energy inside the PWN 
■Bint (tagc) ~ 2.1 x 10 48 erg for our model is similar to theirs 
£int(tagc) = (1 - 2) x 10 48 erg, i.e., PSR J1124-5916 has re- 
leased similar amount of its rotational energy for both mod- 
els. We consider tha t our result is consistent with C05 and 
IWinker et all (|2009l ). Lastly, it should be noted that we do 
not need the low-energy component of the injection spec- 
trum Qi n j similar to the case of GO. 9+0.1 studied in TT11 
(radio-emitting particles are created by adiabatic cooling of 



particles of 7 ^ 7b) 



4.3 Gil. 2-0. 3 



Gll. 2 -0.3 is a composit e SNR at d = 5kpc (|Becker et al.l 
Il985l ; I Green et al.l Il988j ) , and the associated pulsar (PSR 
J1811-1925) has L in = 6.8 x 10 M erg s" 1 and r c = 
23kyr (n = 3) (|Torii et al.l Il999l ; iGavriil et akl |2004| ). 
A mo del spectrum wi t h obs ervations is shown in Fig- 
ure G3 |Xh aroman et all (|2007T) found a TeV 7-ray source 
(HESS J1809-193) in the vicinity of Gll.2-0.3. Although 
some discussions about an association between them have 
been made, their a s sociat io n is consid e red u nlikely by 
iKargaltsev fc Pavlovl (|2007l ); iDean et al.l (|2008l ). We con- 
clude that they are not associated from Figure [5] We also 
confirm that the hard X-ray flu x at ~ 10 19 Hz is not ex- 
plained by the PWN emission as IDean et"al] (|2008l ). They 
concluded that it is a contribution from PSR J1811-1925. 

No other broadband spectral study of Gll.2-0.3 is 
found in the literature to our knowledge. Despite the large 
t c of PSR J1811-1925, Gll.2-0.3 is thought as a young 
SNR/PWN system because of its almost spherical morphol- 
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Figure 5. Model spectrum of Gil. 2-0. 3. The observational data, 
the upper limits and the sensitivity of CTA ( 50 hours) a r e also 
plotted. The observed data are take n from iTam et al, (2002) 
(radio). iRoberts et al.l ||2003)) (X-ray). iDean et al.H20oj) (X-ray 
from PSR J1811-1925) and Uharonian et al.l i2007l ) (7-ray from 
HESS J1809-193). Used parameters are tabulated in TableQ] 
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Figure 6. Model spectrum B0540-69.3. The observational data, 
the upper limits and the sensitivity of CTA ( 50 hours) are also 
plotte d. The observed data are tak e n from [M aiichcstc r~et al.l 
119931 ) (radio). iGallant fc Tuffs! lll998l);ISerafimovich et al.l (120041) 
(infra red and optical), Hiravama et al. I | |2002|) : I Kaaret et al,l 
|200ll ) (X-ray), ICampana et all fcOOsF (combined flux of PSR 
and PWN), and iKomin et al.l ll2012T) (-y-rav from N157B). Used 
parameters are tabulated in Table [T] 



ogy of the PWN and a small off-set of the pulsar from 
the center of the SNR. Moreover, Gil. 2-0. 3 is suspected 
to be associat ed with a historical SN in A.D. 386, i.e., 
iag C ~ 1.6kyr |Clark fc Stephenson! 1977) and obse rvations 
of expan sion rates in radio ijTam fc Roberts! 12003! ) and in 
infrared jKoo et al.ll2007T I support their association. C05 ob- 
tained £ a 



1.6kyr. Our result i agc = 2.0kyr is consistent 



with these age estimates. 



4.4 SNR B0540-69.3 



B0540-69.3 is a composite SNR in the Large Magellanic 
Cloud (LMC, d ~ 50kpc), and the associated pulsar (PSR 
B0540-69) has L spin = 1.48 x 10 38 erg s" 1 , r c = 1.7kyr and 



n ~ 2.1 jCradari et al.ll201ll ; [Livingstone et al.ll2005l ). Be- 
cause the ISRF in LMC is uncertain, we tentatively use val- 
ues around the Sun (Um, (7opt) = (0.5eV, 0.5eV). A model 
spectrum with obse rvations is shown in Figure [6] Although 
IKomin et ail (|2012l ) detected TeV 7-rays fr om N157B close 
to B05 40-69.3, B0540-69.3 was not detected. ICampana et"al] 
([2003) showed a combined flux of the PWN and pulsar 
in hard X-rays and argued that the PWN emission con- 
tributes about 75 — 8 %. Our result seems consistent with 
ICampana et al.l l|2008l ). Although we fail to fit infrared and 
optical observations, these fluxes seem to have consider- 
able uncertainties because a spatially varying background 
makes background s ubtraction uncertain as discussed by 
IWilliams et all (|2008l ). 

No other broadband spectral study of B0540-69.3 is 
found in the literature to our knowledge. C05 obtained 
ta, S e ~ 0.8k y r wh ich is consistent with our calculation. 



iMorse et all (120061 ) found a line-width of filaments 
3300km sec" 1 (corresponding expansion speed is a half of 
this value) which is also consistent with our calculation. 
However, they conclude t asc ~ 1.2kyr because faint mate- 
rials extend (~ 8" ) slightly larger than the PWN ~ 5" . We 
consider that i agc ~ 0.8kyr in our model is acceptable. If 
the ISRF around B0540-69.3 is close to or larger than the 
assumed values (Um,U PT) = (0.5eV, 0.5eV), B0540-69.3 
may be detected by future 7-ray observations. 



5 DISCUSSIONS 

In Sections [3] and [4] we found that all the young PWNe 
can be fitted by a s imilar value of the fracti on parameter 
77 ~ a few x 10~ 3 . iBucciantini et al] (|201ll ) also studied 
spectral evolution of young PWNe, but their parameter cor- 
responding to rj in our model takes a broad range of values. 
This is because a magnetic evolution in their model is differ- 
ent from ours (Equation (|A3|) ). Both their and our models 
are not rigorous treatments based on magnetohydrodynamic 
equations, and then we do not discuss which is a better ap- 
proximation. However, because we presume it plausible that 
properties of all young PWNe are similar to each other, in 
the discussion below, we do not consider the case when 77 
significantly differs from a few x 10 -3 , i.e., the following 
discussions are based on the parameters in Table [T] We es- 
pecially focus on the differences and similarities of TeV and 
non-TeV PWNe. 

We first discuss the detectability of TeV 7-rays from 
young PWNe. We search for correlations between all the 
parameters and the detectability in TeV 7-rays. Figure [7] 
shows the correlation between the total injected energy E to t 
versus the ISRF energy density in infrared ?7ir. TeV and 
non-TeV PWNe are shown by circles and triangles, respec- 
tively. As seen from Figure [7] TeV PWNe are located on the 
top right region (large values of i?tot and Um) and non-TeV 
PWNe are located on the bottom left region (small values 
of Etot and Um,). This behavior indicates the importance of 
total particle energy stored in PWNe and the energy den- 
sity of target photons of inverse Compton scattering for the 
de tectability of T eV 7-rays from young PWNe. As claimed 
bv lHelfandl (|2007f) . the local ISRF energy density is impor- 
tant for the detectability of TeV 7-ray. The Crab Nebula 
(the lowermost circle in Figure [7} is in the lowest ISRF en- 
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Figure 7. The correlation between the total injected energy -Etot 
versus the ISRF energy density in infrared Cir, of young PWNe. 
TeV PWNe are shown by filled circles except for Kes 75 (open 
circle) and non-TeV PWNe are by filled triangles. Young TeV 
PWNe arc located on the top-right region, and young non-TeV 
PWNe are located on the bottom-left region. 
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Figure 8. The correlation between the initial rotational energy 
Loto versus the magnetic energy B^R'l/Q of the central pulsars. 
TeV PWNe are filled circle except for Kcs 75 (open circle) and 
non-TeV PWNe arc filled triangle. 



vironment, although its 7-ray emission is dominated by SSC 
component. Gil. 2-0. 3 (the uppermost triangle in Figure [7]) 
is TeV-undetected although it is located in a high ISRF en- 
vironment. The positions of these two objects in Figure 7 
suggest that the local ISRF alone does not determine the 
TeV 7-ray detectability. For other combinations of the pa- 
rameters, such as the distance d (determining the observed 
flux), the radius Rpwn (determining the degree of adiabati c 
cooling), or r c and L sp - ln (studied bv lMattana et al]|2009t ). 
we do not find any significant correlations, and we do not 
show them here. 

The properties of central pulsars may be different from 
TeV to non-TeV young PWNe. Following TT11, we charac- 
terize each pulsar with the initial rotational energy and the 
magnetic energy. In Figure [8j we plot the correlation be- 
tween the initial rotational energy Lqto versus the magnetic 
energy of the pulsar E B = B»H»/6, where B, <x p 1 / 2 p 1 / 2 
(assuming magnetic dipole radiation) and R* = 10 6 cm are 



the surface dipole magnetic field and the radius of the pul- 
sar, respectively. We find no significant difference between 
the central pulsars of young TeV and non-TeV PWNe. Ex- 
cept for Kes 75, the initial rotation and the magnetic en- 
ergies are uniformly distributed between Lqtq ~ 10 48_50 erg 
and E B ~ 10 42_44 erg (2 x 10 12 - 3 x 10 13 G for the surface 
dipole magnetic field strength). 

For the parameters of a particle distribution at injec- 
tion Qmj(7,£), while we found that P2 = 2.5 is common 
to TeV PWNe, we find no typical value for each parame- 
ter anymore when we include non-TeV PWNe. It may be 
because that spatially varying X-ray spectral indices which 
are observed from most PWNe do not guarantee a fitting 
by a one-zone model. It is not clear why TeV PWNe can be 
fitted by a common value of P2 = 2.5. On the other hand, 
radio and optical emissions are more spread in entire PWNe 
than X-ray, i.e., showing approximately one-zone behaviors, 
and then the varieties of pi and 7b seems real. Interestingly, 
G292.0+1.8 does not need the low-energy power-law com- 
ponent at injection as TeV PWN GO. 9+0.1 also shows the 
same characteristic (TT11), but we do need the low-energy 
component for other young PWNe. Lastly, following TT11, 
we can estimate a lower limit of the multiplicity k and a 
upper limit of the bulk Lorentz factor P w of a pulsar wind 
(listed in Table [1]). These parameters are also common to 
TeV and non-TeV PWNe. 



6 CONCLUSIONS 

We applied the spectral evolution model developed by TTU 
to young non-TeV PWNe, 3C58, G310.6-1.6, G292.0+1.8, 
Gil. 2-0. 3 and B0540-69.3 to study differences between non- 
TeV and TeV PWNe. Observed TeV 7-ray flux upper limits 
of these non-TeV PWNe give us only upper or lower limits 
of the model parameters. However, we can estimate plau- 
sible values of parameters from other observational or the- 
oretical studies which are independent off our calculations. 
For the case of 3C58, the fitting parameters are significantly 
constrained and a lower limit of the fraction parameter is 
77 ^ 3.0 x 10 -3 with t ag c = 2.5kyr. Combined with the age 
estimated from the expansion velocity of filamentary struc- 
tures, we consider that r\ ~ 3.0 x 10 -3 with i agc ~ 2.5kyr 
is reasonable for 3C58. This value of fraction parameter is 
similar to that of young TeV PWNe (rj ~ a few x 10~ 3 ). For 
other young non-TeV PWNe, we also find that the fraction 
parameters of r\ ~ a few x 10~ 3 is consistent with observa- 
tions although a wide range of values is nominally allowed. 
Our model alone weakly constrains the fraction parameters 
with a smaller value of a lower limit. The comparison be- 
tween other observational and theoretical studies makes the 
fraction parameters of young PWNe being r\ ~ a few x 10 -3 
acceptable. 

For the detectability of young PWNe in TeV 7-rays, we 
find the detectability is correlated with both the total in- 
jected energy E to t and the ISRF energy density in infrared 
(7ir . This behavior indicates importance of the total parti- 
cle energy stored in PWNe and the energy density of target 
photons of inverse Compton scattering for the detectability 
of TeV 7-rays from young PWNe. We do not find any signif- 
icant correlations for other combinations of the parameters. 

The central pulsar properties, i.e., the initial rotational 
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energy Lqto and the magnetic energy Eb of the pulsar, are 
not significantly different between TeV and non-TeV PWNe. 
The values are distributed in a range of Lqtq ~ 10 48 _50 erg 
and E B ~ 10 42 ~ 44 erg. 

The fitted break energy 7b and the high energy power- 
law index at injection P2 are widely distributed, 7b ~ 10 4-7 
and p 2 ~ 2.5 - 3.0. Interestingly, G292.0+1.8 does not 
need the low-energy component at injection to reproduce 
the observed radio spectrum. This behavior is the same as 
GO. 9+0.1 studied in the previous paper, i.e., the adiabatic 
cooling of the high-energy component reproduces radio emis- 
sions from the object. 

The radio flux evolution of 3C58 calculated in our model 
is about — 0.04%yr _1 which is four times smaller than that 
of the Crab Nebula. 3C58 is the only one PWN whose radio 
flux evolution are obtained, except for the Crab Nebula. The 
latest WMAP observation shows that the flux evolution is 
almost zero within error and then we conclude that the radio 
flux evolution calculated in our model is consistent with the 
WMAP observation. 
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The particle injection Qinj(7, t) is assumed to follow a broken 
power-law characterized by five parameters 7 m i n , 7b, 7max, 
pi and p2, which are the minimum, break and maximum 
Lorentz factors and the power-law indices at the low and 
high energy ranges of the injection spectra, respectively. The 
particle injection Qmj(7,i) has a form of 



>inj(7,i) 



Qo(t)(7/7b) _pl for 7min S 7 ^ 7b , 
Qo(i)(7/7b)" P2 for 7 b S 7 S 7u 



(A2) 

where 7 is the Lorentz factor of relativistic particles. We 
require that the normalization Qo{t) satisfies (1 — rj)L(t) = 
j7max Q in j(-y ; t)ym e c 2 d'y, where m e and c are the mass of an 
electron (or positron) and the speed of light, respectively. 
The cooling term 7(7, t) includes the radiative coolings (syn- 
chrotron radiation and IC/ISRF) and the adiabatic cooling. 
Lastly, we need to specify the magnetic field evolution, which 
is assumed to be given by magnetic energy conservation, 



y(flpWN(t)) 3 



(B(t)f 
8n 



[ r)L(t')dt' 
Jo 

r]E tot (t), 



(A3) 



where E to t(t) is the integrated spin-down power at a time t. 
Although Equation ()A3[) neglects adiabatic losses, the mag- 
netic field also suffers from winding and stretching process. 
Uncertainties of one-zone models may be larger than these 
corrections. Note that some j ustifications of Equat ion (1 A3f) 
are discussed in Section 2.2 of iTanaka fc Takaharal (|2010l ). 



APPENDIX A: SPECTRAL EVOLUTION 
MODEL OF YOUNG PWNE 

We have modeled spectral evolution of young PWNe. The 
radiation processes are synchrotron radiation and inverse 
Compton scattering off synchrotron radiation (SSC) and 
off the ISRF (IC/ISRF). The ISRF has three components 
(the cosmic microwave background radiation (CMB), in- 
frared photons from dust grains and optical photons from 
stars) . The infrared and optical components of the ISRF are 
parametrized by the energy densities Um and Uopt which 
vary with the distance from the Galactic center r and the 
height above the Galactic plane z. 

We consider a PWN which is a uniform sphere expand- 
ing at a constant velocity i>pwn (the radius of the PWN 
fipwN(t) = DpwNi). The PWN contains a magnetic field 
and a population of accelerated electron-positron plasma, 
which are injected from its central pulsar. The spin-down 
power L sp i n (f) of the central pulsar is divided into the mag- 
netic power ^L sp i n (t), and the particle power (1 — r])L sp i n (t), 
where r\ (0 ^ 77 ^ 1) is the fraction parameter. Evolu- 
tion of the spin-down power is given by L ap i n (t) — Lq[1 + 
(t/To)]~^" +1 ' / '' n_1 ' , where to is the spin-down time and Loto 
is the initial rotational energy. We need four values to specify 
Lspin(t), the current pulsar period P, its time derivative P, 
braking index n and the current age i age . The characteristic 
age r c = P/2P is described as r c = (n — l)/2 • (to + £ ag e). 

For evolution of a particle distribution A r (7,f), we use 
the continuity equation in the energy space, 

|W(7, t) + .g- (7(7, t)N{j, t)) = Q inj (7, t). (Al) 
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